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The magnetization detection and switching of an ultrasmall Stoner nanograin in a non-local spin 
valve (NLSV) device is studied theoretically. With the help of the rate equations, a unified descrip- 
tion can be presented on the same footing for the NLSV signal that reads out the magnetization, 
and for the switching process. The setup can be viewed as that the grain is connected to two 
non-magnetic leads via sequential tunneling. In one lead, the chemical potentials for spin-up and 
-down electrons are split due to the spin injection in the NLSV. This splitting (or the spin bias) is 
crucial to the NLSV signal and the critical condition to the magnetization switching. By using the 
standard spin diffusion equation and parameters from recent NLSV device, the magnitude of the 
spin bias is estimated, and found large enough to drive the magnetization switching of the cobalt 
nanograin reported in earlier experiments. A microscopic interpretation of NLSV signal in the se- 
quential tunneling regime is thereby raised, which show properties due to the ultrasmall size of the 
grain. The dynamics at the reversal point shows that there may be a spin-polarized current instead 
of the anticipated pure spin current fiowing during the reversal due to the electron accumulation in 
the floating lead used for the readout of NLSV signal. 

PACS numbers: 73.23.-b, 85.75.-d, 72.25.Hg, 85.35.-p 



I. INTRODUCTION 

Current-induced magnetization reversal had attracted 
considerable interest due to its fundamental significance 
in understanding interplay between magnetism and elec- 
tricity as well as potential applications in magnetic 
memories i^i^i^i^i^ii As the scale of the ferromagnetic 
nanograins goes down to only several nanometers j^iiSiiii 
many theoreticali^ii^ii^ii^iiiii^ii^i^aiSi and experimen- 
tal works22i22ii^i2^ were inspired to address the current- 
induced magnetization reversal in these small structures. 
By far, most studied setups were multi-layer or nanopillar 
structures with vertical geometries, in which spins are al- 
ways carried along the flowing of charge current. Usually, 
the critical current density as high as 10^ — 10^ A/cm^ 
is required to induce a reversali^ Considering such high 
density of current flowing through each nanograin, when 
a huge amount of nanograins are integrated in large scale, 
spurious effects such as Joule heat, current-induced mag- 
netic field, and noise are not ignorable. 

A possible solution is to use pure spin current, in 
which the same amount of spin-up and -down currents 
flow along opposite directions, yielding no net electric 
current. By far, one of the most promising designs to 
realize considerable pure spin current is the non-local 
spin valve (NLSV) devices21.^i29^i2i.32^|34 ^^^^ j^^^g^^i 

geometry^2^ Recently, a NLSV was reported to reversibly 
switch the magnetization of a ferromagnetic particlei^i^ 
As shown in Fig. [ija), a typical NLSV includes a big- 
ger fixed and a smaller free ferromagnets (denoted by 
the shadow areas) embedded on the left and right sides, 
usually referred as injector and detector, respectively. 
By driving a current Ic through regions (2), (1), and 
(3), spins can be injected from injector (1) to produce 
a nonequilibrium spin accumulation in nonmagnetic re- 



gion (4). This spin accumulation exhibits as a splitting of 
chemical potentials for spin-up and -down electrons (spin 
bias or spin voltage), as shown in Figs. [3ja) and [3Kb). 
While the spin-polarized current flows only in the loop 
formed by regions (2), (1) and (3), the nonequilibrium 
spin accumulation in region (4) diffuses to the right ac- 
companied by a pure spin current flowing toward the de- 
tector. In this way, net charge current is prevented from 
flowing directly into the grain. In the response of Ic, the 
magnetization of the detector can be read out by mea- 
suring the voltage difference AV between (5) and (6), re- 
ferred to as the NLSV voltage, which is usually estimated 
by the spin diffusion equationi^Sii^Si^i^i^ By applying Ic 
exceeding a critical value, the magnetization of the detec- 
tor can be switched reversiblyi^SiSl This process is usually 
described by Landau-Lifshitz-Gilbert (LLC) equation.''— 

By far, the detectors are fllms of size of 100 - 1000 
nm. Considering commercial charge devices already work 
well at these scales, the ultimate goal of utilizing pure 
spin degree of freedom is to replace the charge devices 
at only several nanometers. Besides, because the LLC 
equation pre-assume the magnitude and polarization of 
currents that flow through the grain as input parameters, 
the counteraction of the detector on the currents is not 
taken into account. As a result, important information 
could be missing, e.g., whether the pure spin current is 
still a pure spin current after flowing into the detector. 

In this work, we study a NLSV device, in which the 
usual film detector is replaced with one or multiple well- 
separated cobalt nanograins embedded in insulator, as 
those in Ref. dj. The nanograin is much smaller in size 
(^ 10'^ atoms) than the films and at much lower tempera- 
tures {'^ 20 mK). Such small nanograin can be viewed as 
a Stoner particle whose fcrromagnetism comes from the 
exchange interactions between itinerant electrons inside 
it, and thereby can be manipulated by exchanging an- 
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gular momenta with the electrons that tunnel through 
itJ^iii The grain is modeled as coupled to two non- 
magnetic leads via quantum tunneling. In one lead, the 
chemical potentials for the spin-up and -down electrons 
are split, due to the spin injection in the NLSV. The 
small size of the nanograin allows us to model the NLSV 
signal and the magnetization switching within one set of 
rate equations. Besides, the previous knowledge of the 
Co nanograin from experimenta^i^ i^^'^^ and theories^i^ 
allows us to perform a realistic evaluation. 

This work focuses on two aspects: (i) the possibility 
of the detection and switching is evaluated using the pa- 
rameters extracted from the previous NLSV"^®'"^^ and Co 
grain^iiSiiii experiments, (ii) Such small grain is sub- 
jected to strong Coulomb and magnetic blockades,^® how 
these blockades determine the critical conditions for the 
reversal, e.g., the critical driving current Ic, the gate volt- 
age V^, and the spin bias Vs- 

We find that: (i) The numerical evaluations using re- 
alistic parameters from the recent NLSV— i2i and the 
cobalt grain experiment o^^'^°i^^ show that it is possible to 
employ the NLSV device to detect and switch the magne- 
tization of a ferromagnetic nanograin under the present 
experimental conditions, (ii) Under /c, the NLSV signal 
can also be detected in the sequential tunneling regime to 
read out the magnetization of the grain, and interpreted 
from a microscopic view. Interestingly, if the majority 
band of the grain is favored to participate in the electron 
transport, the sign of the NLSV signal turns out to be 
just opposite to that if the the minority band is preferred 
to conduct electrons. In the presence of an angle 9 be- 
tween the easy axis of the grain and the spin-quantization 
direction of the lead, the NLSV signal is proportional to 
cos 9 and vanishes when — 7r/2. (iii) Under Ic ex- 
ceeding a critical value, the magnetization of the grain 
can be switched reversibly. The critical current Ic re- 
quired for the magnetization switching is determined by 
the gate voltage Vg and the spin bias Vs at which both 
the Coulomb and magnetic blockades in the grain are 
lifted. By choosing suitable gate voltage, the critical Vs 
needed for switching can be minimized to 2Kj^, where 
Kn is the volume-independent anisotropy of the grain. 
Besides, the transient current during the magnetization 
reversal may be a spin-polarized current instead of the 
anticipated pure spin current, due to the electron accu- 
mulation or drainage in the floating lead used for the 
NLSV measurement. A possible solution is to remove 
the floating lead. 

The paper is organized as follows. In Sec. [TTl the 
model and theoretical formalisms are introduced. In Sec. 
IIII[ the microscopic NLSV signal in the sequential tun- 
neling regime is described in detail. In Sec. llVi the 
magnetization switching under large injection current Ic 
is presented. Finally, a summary is given in Sec. |Vl 



II. MODEL AND THEORETICAL 
FORMALISMS 
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FIG. 1: (a) The schematic of our non-local spin valve (NLSV) 
device. It can be viewed as a combination of the setups em- 
ployed in Ref. and Ref. @. (b) The device is separated into 
7 regions. (1) is a fixed ferromagnetic injector of thickness d. 
The ferromagnetic grain on the right is to be manipulated 
and detected. The distance between (1) and the grain is L. 
" ^ " in each region defines the local origin and positive 
direction of coordinate, (c) The grain and its nearby regions 
is modeled as a Stoner particle coupled to two non-magnetic 
lead via tunneling through two barriers. Driven by the spin- 
polarized current Ic, spins injected from (1) accumulate and 
induce a splitting (Vs) of the spin-up and down chemical po- 
tentials [/i''|,''^(I/)] at the (4)/(5)/grain interface. Vg is the 
gate voltage applied to the insulator [see Fig. 1(a)] surround- 
ing the grain, much like the technique in Ref. [y. In this way, 
it can be used to tune the energy levels inside the grain with 
respect to chemical potentials of (4), (5), and (6) capacitively, 
while without inducing direct current between the gate and 
the grain. AV^ measures the voltage between (5) and (6), and 
is defined as the NLSV voltage, /i'^' (cxa) is the chemical po- 
tential at the voltmeter side of regions (5), and is set as the 
energy zero point throughout the paper. 
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A. General survey of our setup 

The device we study is shown in Fig. [Ija), which can 
be divided into 7 regions, "(l)-(6)" and "grain", as shown 
in Fig. [Hb). The different regions of the device are 
separated into three parts, and modeled by different for- 
malisms, depending on their sizes and positions: 

(i) the first part consists of regions (l)-(5). Their sizes 
are comparable to their spin diffusion lengthes, thus the 
spin transport in these regions are governed by the spin 
diffusion equationsi^i^'^i^i^ The magnetization of in- 
jector (1) is assumed to be fixed. 

(ii) the second part is the smaller ( approximately 
several nm) grain, which is described as a Stoner 
particler^i^ coupled via sequential tunneling to two non- 
magnetic leads. For convenience, we call them lead (4,5) 
and lead (6), respectively. Lead (4,5) and lead (6) are 
defined as where regions (4) and (5) and region (6) con- 
nect the grain, respectively. We assume that there is no 
direct tunneling between lead (4,5) and lead (6), i.e., the 
only possible connection between them is via the grain. 

(iii) The third part is floating region (6). Because of 
the voltmeter, it is an open circuit between regions (4) 
and (5) and region (6), i.e., at steady state, there is no 
current flowing from regions (4) and (5) across the grain 
to (6), due to a voltage difference between (6) and (5). 
Throughout the paper, we define the chemical potential 
/i^^^ (oo) at the voltmeter side of region (5) as the energy 
zero point. The chemical potential of region (6) is de- 
noted as /i*^®^ . The difference between and (oo) is 
denoted as AV = /z^^^ - (oo), where AV is the NLSV 
voltage. AV will be determined through self-consistent 
calculation. 

In a word, our model can be viewed as a combina- 
tion of the setup used in two kinds of experiments, i.e., 
the NLSV'^^ and the transport through ferromagnetic 
nanograins.^'^ Besides: 

(i) Different from usual ferromagnetic nanograins 
setupsj^i^ the chemical potentials for spin-up and -down 
electrons fj,'"y^{L) are split at the (4)/(5)/grain interface, 
as shown in Fig. [IKc). The splitting is called the spin 
bias, and denoted as This spin bias is 

induced by the spins injected from region (1) to region 
(4). Later in Sees. IIIII and llV Al we will see that Vs is 
crucial to the detection of NLSV signal and the magneti- 
zation reversal, thus we will first evaluate its magnitude 
in the following subsection. 

(u) DiflFerent for usual NLSVSfS^i^ The ultrasmall size 
of the grain allows the NLSV signal, the magnetization 
reversal dynamics, and the interaction between the cur- 
rents and the grain to be described within one set of rate 
equations Eq. 

For a better understanding of our setup, the calcula- 
tion steps of the magnetization and the NLSV signal as 
functions of Ic is presented in Fig. O This section is 
organized as follows. In Sec. Ill Bl Vs will be estimated 
with realistic experimental parameters. In Sec. Ill C[ the 



description of the grain will be introduced. In Sec. IIIDI 
the transport between the grain and its leads will be de- 
scribed by the rate equations. The details of theoretical 
descriptions can be found in Appendixes EUCl and will 
be specified in the following subsections. 
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FIG. 2: The calculation steps of the magnetization of the 
grain and the non-local spin valve signal -Rnlsv as functions 
of 7c, the driving current. Rectangles represent calculation 
processes. Parallelograms represent inputs or outputs of the 
calculations. The region (1), (2), (3), (4) and (5) in Fig. [1] 
are described by the spin diffusion equation given by Eq. 
The grain and its nearby regions are described by the rate 
equations given by Eq. (|12|) . The diamond represents the 
self-consistent calculation of /i''^', the chemical potential of 
region (6). 



B. Spin diffusion transport in regions (l)-(5) and 
estimate of spin bias Vs 



In this section, we will estimate the magnitude of 14, 
which is defined as the splitting of spin-up and -down 
chemical potentials at the (4)/(5)/grain interface, 



y, ^ i,f{L) - ^(L) ^ ^^^^(0) - Mr(0). (1) 

Following the previous literatures, the spin and charge 
transports in regions (l)-(5) are described by the stan- 



,(4) 



(5), 



,(5), 
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dard diffusion equation^i^iiaiilii^ 



V (cti^i +crxAii) = 0, 
V'(MT-m) = ^ 



Mi), 



(2) 



where A is a phenomenological spin-diffusion length that 
can be measured experimentally.^^ Usually, the spin- 
diffusion length in normal metal is much longer than that 
in ferromagnetic metal, i.e., Xn ^ Xp. In this work, 
we use Xn = 1000 nm, and Xp = 5 nm.^^ is the 
conductivity for a electrons. For ferromagnetic metal, 
(Tf = ap{l + a)/2 and fT| = ap{l — a)/2, where ap is 
the total conductivity of the ferromagnetic metal, a is 
the bulk polarization defined as a = (cr| — o'i)/(o'T + "'i)- 
For normal metal, cri = (T| = ctn/'^- In this work, we 
use (TAT = 8.8 X 10'^n-^m-\ ap = 9.8 x lO^fl-^m-^, and 
a = 0.2,^ 

We calculate the chemical potentials in regions (l)-(5) 
of Fig. HUb), where regions (2), (3), and (5) can be re- 
garded as semi-infinite. By assuming that the same cross- 
section area Aq perpendicular to the current in each re- 
gion, the solutions to Eq. ^ can be simplified to be 
one dimensional. Following Jedema et al.,'^^ the general 
solutions to the 5 regions as functions of position x are 
given by 
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where the origin and positive direction in each region are 
locally defined for a concise form of the general solutions, 
and indicated by " — > " in Fig. [Tfb). jc = Ic/A^, where 
Ic is the spin-polarized total current flowing through re- 
gions (2), (1), and (3). e is the electron charge. The 
9 coefficients A, B, C, D, E, F, G, H, and K will be 
determined by boundary conditions. 

The boundary conditions are given as follows: 
(i) The current polarization usually loses when flowing 
through an interface from the ferromagnetic to the nor- 
mal side, due to, e.g., the spin-dependent scattering. We 
take this loss into account by phenomenologically intro- 
ducing an efficiency parameter 7 G [0, 1], 
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(4) 



where — —{ua / e)dx^.(j, is the current density for the 
spin-cr electrons. By combining Eq. ([4]) and the conser- 
vation of the total current 



(5) 



the boundary conditions at the ferromagnetic/normal in- 
terface are then given by. 



7 + 1 . 7 — 1 
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(6) 



(ii) The chemical potentials of each spin components 
are continuous at each interface. For simplicity, we do not 
explicitly include the spin-dependent chemical potential 
drops caused by the interface resistance. Its destructive 
effects, in particular on the reduction in the spin bias, will 
be approximately accounted by considering a relatively 
small injection efficiency 7. 

(iii) At the (4)/(5)/grain interface, the spin-up and 
-down currents flowing into the grain are l'^ and — I*^, 
respectively, i.e., (how this boundary condition is derived 
can be found in Appendix [X]) 



Au^-^dJ^)^{L)^±eI^ 



(7) 



where we assume a pure spin current flowing from (4) 
and (5) into the grain, which however may not be true 
as we will see in Sec. IIVBI However, this deviation is 
neglected because I'^ is too small to affect K , as we will 
see in Sec. IIVCI Actually, because J*^ is ignorably small, 
we simply neglect it in our numerical calculations, though 
we include it in equations explicitly. 

With these boundary conditions, the coefficients from 
A through K in Eq. ^ are readily found. Figs. [H^a) 
and EJb) show the spin-resolved chemical potentials /i| 
and /i| in regions (l)-(5). For a clear demonstration of 
the splitting between spin-up and -down chemical poten- 
tials, we choose a — 0.9 only in these two figures. For 
the realistic evaluations in the rest part of this work, we 
choose a = 0.2, as estimated by the experimentsi '^^'^^ All 
the conductivities, spin-diffusion lengthes, and the bulk 
polarization used for the evaluation are extracted from 
the experiment AaXa.^^ and given in Fig. [31 Driven by 
Ic, the spins injected from injector (1) induce a splitting 
between ^| and at interfaces (l)/(2) and (l)/(3)/(4). 
Because injector (1) is sandwiched between two normal 
metals, ^| and //j^ in (1) cross with each other at the 
center and split oppositely on the opposite sides. In each 
region of (1), (2), and (3), because of the charge current 
Ici A^T ^^'^ Mi not only split but also demonstrate steep 
slopes with the same trend. Because of no net charge cur- 
rent, fjL'\ and 111 in region (4) just diffuse to the right side 
in opposite gradients (not obvious here because L <^ Xjq)- 

By using Eq. ([1]), the analytic expression for Vs is 
found out as 
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Xp a {Ml 
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where 
Ml 



2cosh(-^) + 2Afosinh(-^), 



M2 = 6A/oCosh(-^) + 2(2M2 + l)sinh(-^^), 

Xp Ap 



Mq = 



ajyXp 



(9) 



Figures. [3lc)-(e) show the spin bias Vs as a function of Ic 
and 7. Later we will see in Sees. IIV Bl and IIV CI that 
is ignorably small, thus hardly changing Vs- Thus, we let 
= 0. As we see, Vs can be as large as ~ 0.1 meV for 
the present device and parameters. Vs increases linearly 
with /c, while logarithmically with 7. Vs changes sign 
with the current L- 
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FIG. 3: (a)(b): The chemical potentials /X| (marked by |) 
and fii (marlted by J,) in regions (l)-(4). They change with 
the direction of Ic, the spin-polarized current flowing through 
regions (1), (2), and (3). Vs is defined as jif — fii at the 
(4)/ (5) /grain interface. The horizontal dashed line indicates 
/i'^'(cx)), the energy zero point throughout the paper, (c) Vg 
vs Ic for different values of spin injection efficiency 7. (d) 14 
vs 7 for different Ic. (e) Vs (in unit of meV) as a function of 
Ic and 7. The parameters for (c)(d)(e): d = 20 nm, L = 270 
nm, Aq = 170 x 65 nm^, Af = 5 nm, Ajv — 1000 nm, ap = 
9.8 X lO'' n-^m-\ aN = 8.8 x 10^ fi^^m^S a = 0.2, /° = 0. 

In the following, for a conservative and realistic sim- 
ulation, we will always choose a relatively low injection 
efficiency 7 = 0.1 and assume = 0. Since the above 



boundary conditions only consider the conservation of 
current density, our results are valid when assuming the 
identical cross-section area in each part of the setup. One 
should consider different cross section area in different re- 
gions and current conservation for a more general case4^ 



C. IVIany-body states of the ferromagnetic 
nanograin 

In this work, we will describe the ferromagnetic 
nanograin by using the minimal possible model^i^ 
proposed to describe the experiment transport spec- 
tra through cobalt nanograiui^ This model was 
also adopted to discuss spin-polarized current-induced 
relaxation and spin torque in the ferromagnetic 

16 1 7 44 

nanogramj °i 

A full and detailed description of this model can be 
found in Appendix |B] Simply speaking, at low temper- 
atures, the grain can be described by the many-body 
states \N,S,Sz), where N is the total electron number 
inside the grain, S and Sz are the magnitude and the 
z-component of the total angular momentum S of the 
grain, respectively. 

In the following, we will focus on two branches of 
states. The first branch is 



\No, So, Sz), 



(10) 



i.e., there are N — Nq electrons inside the grain, and the 
magnitude of the total angular momentum of the grain is 
S = So- Besides, since Sz G [—5*0, So], there are 25o + 1 
states in this branch. The second branch is obtained 
by adding an extra electron to the minority band of the 
grain with respect to the first branch, so that the total 
electron number increases by 1 and the magnitude of the 
total angular momentum decreases by 1/2. This branch 
is denoted as 



\No + I, So- 1/2, Sz), 



(11) 



where G [— 5o + 1/2, 5o — 1/2], i.e., there are 250 states 
in this branch. 

We refer regions (4)- (6) as the two "leads" connected 
to the grain, one is from regions (4) and (5), the other 
is from region (6). By tuning the gate voltage Vg and 
applying the spin bias Vs (induced by Jc), the energies of 
the two branches presented in Eqs. ([TO]) and pT|) can be 
set to be nearly degenerate with respect to the chemical 
potentials /i|-'*j(L), /i|^|(0), and ii^^K In this situation, 
the electrons in lead (4,5) and lead (6) can be exchanged 
with the grain. Then we can use the polarization of the 
exchanged electrons to detect and manipulate the mag- 
netization of the ferromagnetic grain. 



D. Rate equations in the presence of spin bias 

The evolution of the many-body states of the grain 
by exchanging electrons with the weakly coupled lead 
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FIG. 4: The easy-axis of the grain is set as z direction. 
The grain electrons are quantized along z direction. 

The majority and minority electrons orient parallel and anti- 
parallel with S, respectively. S is the total angular momen- 
tum of the grain formed by all the electrons inside the grain. 
Q is the angle between the easy axis of grain and the spin- 
quantization direction of lead (4,5). 



Qa in Eq. p4| is the angle between the easy axis of 
grain and the spin-quantization direction of lead a. For 
simplicity, we set the easy-axis of the grain as z-axis and 
assume that 0(g) — and ^(4^5) — 6 ^ [0,7r/2], as shown 

in Fig. m In the following, we denote — when 

In terms of Pi , the magnetization of the grain is given 

by 

M = Y,SiPi, (15) 

and the spin-u current flowing from region a into the 
grain are defined as 

i:.=-eJ2iNi-N,)R¥^P,, (16) 
w 

where Ni and S*' correspond to the N and of the state 
l = \N,S,S,}. 

The validity of the rate equations is discussed in Ap- 
pendix [cl 



(4,5) and lead (6) is described by the Pauli rate 
equations/^i^i^ We only consider the sequential tunnel- 
ing regime. Born approximation and Markoff approxima- 
tion are applied, and Ht is treated by perturbation up to 
the second order The rate equation can be expressed 
in a compact form, 

dtPi^Y.^U'Pi'^ (12) 

where < P/ < 1 are the probability to find the state 
I = \N, S, Sz)- The diagonal and off-diagonal terms of the 
coefficient matrix of the rate equations are, respectively, 

Ri'^i ^'^^"'^h Ru = ~'^Ri'i, (13) 

aai I'^l 

where 

Rt};^ = r„cos2^[Kr|d,TlOlV(i5;-i?/'-M?) 

+\m^\l')\^f{El-El,+^,%)] 

sin^ ^-^[\{l'\d,^\l)\^f{Ei - El, - fil) 
+ \{l\dM'f{Ei^Ei,+^,'^_)],iU) 

and one just replaces t with J, and exchanges -t- and 
— to obtain Note that the Fermi distribution 

f{x) = l/[exp(x/A;Br) -I- 1] is spin-resolved. The param- 
eter Ta — 27r^j, |V'fcQip(5(a; — eka) represents the spin- 
irrelevant coupling between lead a £ {(4, 5), (6)} and the 
grain. For simplicity, we assume that 5) = r(g) = F, 
and r are assumed to be independent of the specific 
single-particle level i. The overlapping {l\dia\l') can be 
found by calculating the Clebsch-Gordan coefficients. 



III. NLSV SIGNAL IN THE TUNNELING 
REGIME 

In this section, we will present the microscopic picture 
of the NLSV signal in the sequential tunneling regime. 
When a small current E is driven in the loop formed 
by regions (l)-(3) of Fig. [IJa), the relative alignment of 
magnetization between the two ferromagnets can be read 
out by measuring the voltage difference between regions 
(6) and (5). 

A. 6 — Q case 

We will discuss first the case when = 0. For = 0, 
we denote = /^T/i' results can be easily gen- 

eralized for 7^ in Sec. HITB] The left and right 
columns of Fig. [5] show the cases when the magneti- 
zation of the grain is anti-parallel (AP) and parallel (P) 
with injector (1), respectively. First, we consider the 
AP case, i.e., Sz of the grain ~ — Sq. By putting the 
experimental measurement current^! =250 /iA into 
Eq. ((8]) and assuming 7 = 0.1 and I'^ — 0, we esti- 
mate that a positive spin bias Vs ~ 1.84 ^eV will be 
induced at the interface where regions (4), (5) and the 
grain connect, so that ^^y^{L) = ^^^'j^(O) = ±Vs/2. It 
is natural to assume that the distance between the volt- 
meter and the (4)/ (5) /grain interface ^ Aa? = 1 /xm, so 
the split chemical potentials for spin-up and -down elec- 
trons will decay to only one chemical potential /i^^-* (00) 
on the voltmeter side of region (5). Because region (5) is 
nonmagnetic, the decays of spin-up and -down chemical 
potentials are symmetric, so it can be anticipated that 
/i(5)(oo) = [^^f\L)+^if{L)y2 = [nf{Q) + ^ifm/2. 
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FIG. 5: Using NLSV signal in the sequential tunneling regime 
to read out the magnetization of the grain. (al)-(a3) and 
(bl)-(b3) show the cases when the magnetization of the grain 
is antiparallel (AP) and parallel (P) with injector (1), re- 
spectively. ^^^\oo) and /i'-®-' are the chemical potentials on 
the voltmeter side of regions (5) and (6), respectively, (al) 
When Sz ~ —5*0, the grain favors spin-up current flowing 
from lead (4,5) to lead (6). Because it is an open circuit be- 
tween lead (4,5) and lead (6), the electrons flowing into lead 
(6) can not go anywhere but build up and raise the chemical 
potential . (a2) Once the chemical potential is raised 
to be aligned with /^^^''(L), the current flowing and electron 
building-up stop. Then one can measure a voltage difference 
between fi^^^oo) and which corresponds to the NLSV 

voltage. (a3) The spin-up and -down currents as functions of 
M^'^^-At'^Hoo). The NLSV vohage AV = n'-^^ - n'-^\oo) corre- 
sponds to the vertical dashed line at which both currents are 
zero. (bl)-(b3) can be understood similarly. The parameters: 
= 0, 7c = 250 /xA, 7 = 0.1, 1° = 0, 5o = 100, T = 20 mK, 
AVg = —Kn, and other parameters are the same as those in 
Fig. [31 



This situation is depicted by fi'"y^{L) and /i'^^ 
[5](al). 



(oo) in Fig. 



In the response of this small spin bias, a small current 
will be generated, flowing through the grain. Specifically, 
spin-up current will be favored for the AP case. This 
can be understood with the help of the Clebsch-Gordan 



coefficient. 



(iVo, 5*0, SMnWo + 1, 5*0 - i 5, + i) 



(ji = So, mi = Sz;j2 = ^, 



■ni2 



So — Sz 
2So + 1 



{No,So,SMl\No + l,So~^,Sz-^) 



(ji = .So, mi = Sz;j2 = -, 
\J^So-l,M^Sz-h 



m2 



So + Sz 
25o + 1 ' 



(17) 



apparently, for 5*2 ^ — ^o, the probability for spin-up 
electrons to tunnel through the grain 



So — i^So) 



2So 



1 



1, 



(18) 



is much larger than the probability for spin-down elec- 
trons 



So + i^So) 
25*0 + 1 



0. 



(19) 



In other words, the spin-down current is magnetic 
blockadedJ^ Note that the spin selection rules remain 
qualitatively unchanged even for small fluctuation of Sz 
around —So, as long as 5*0 S> 1. 

As a result, the favored spin- up electrons will flow 
from region (4) through the grain into region (6), if 

/i|"'''(L) > Remember that there is a voltmeter be- 

tween regions (5) and (6), so the electrons tunneling into 
region (6) can not go anywhere but accumulate and raise 

the chemical potential of region (6) until fj,^^^ 
After this accumulation is accomplished, no more current 
will flow and there is finally a stable voltage difference 
Ay between /i^^^ and /j,*^^^(oo), as shown in Fig. [5](a2). 

Fig. [S] (a3) shows how to numerically determine the 
NLSV voltage AT^. One just scan jjS^^ and calculate the 
tunneling current through the grain. AV is then found 
out as at which /i^^^ —jjS^^ (oo) both the spin-up and -down 
currents vanish, as indicated by the vertical dashed line. 
It turns out that in the sequential tunneling regime and 
when = 0, the magnitude of the NLSV voltage is 



lAV 



(20) 



By using Eq. ([8]) and the parameters given in Fig. [5l for 
Ic = 250 /iA, the NLSV voltage IW^=° is found out to 
be about 0.92 /.iV. The NLSV signal is thereby equal to 



p6»=0 
-n-NLSV 



AF 



0.92 ^V 
250^A 



3.7mf]. 



(21) 
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Similarly, as shown in Figs. [5] (bl)-(b3), the parallel 
(P) case favors that spin-down current flowing from re- 
gion (6) through the grain to regions (4) and (5), also 
due to the same spin selection rules Eq. pT|) . This 
current will drain the electrons in region (6) and lower 
^(6) until fi^^^ ^ fJ.f\L). Therefore, for the P case, 
Ay^=" ^ -0.92/^eV and i?NLsv = -3.7ml7, right op- 
posite to the AP case. 
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FIG. 6: NLSV signal when 61/0. (a) when 61/0 and the 
magnetization of the grain is Sz = —So (AP), the spin-up cur- 
rent can flow simultaneously from lead (4,5) to lead (6) with 
relative probability cos^ |, and from lead (6) to lead (4,5) 
with probability sin^ |. (b) when 9 and the magnetiza- 
tion of the grain is Sz ~ So (P), spin-down current can flow 
simultaneously from lead (4,5) to lead (6) with relative prob- 
ability sin^ |, and from lead (6) to lead (4,5) with probability 
cos^ |. (c) The NLSV signal as a function of 9 for the P (A) 
and AP (V) cases. ±cos6' are also plotted for comparison. 
Other parameters are the same as Fig. \S\ 

When 9^0, the chemical potentials of lead (4,5) now 

are denoted as and ij!"f^{L), respectively. We still 

take Sz ^ —So for example. As shown in Fig. dja), when 
0/0, the favored spin-t current can flow not only from 
fJ^'+'iL) to /.t^^-*, but also from /x^^^ to /x^^(L), according 
to the rate equations Eq. ([H]) . These two currents are 



denoted as / 



and / 



respectively. 



According to the rate coefficients Eq. (|14p . the spin-up 



electrons are related to by cos'^ |, and to p^^L) 

by sin^ |. Besides, lead (6) is not spin-dependent, so 
the polarization and magnitude of the current between 
the grain and (6) just follow those between (4) and the 



,(4), 



grain. Therefore, 



T 2 

/ M> oc cos — , 



rT 



2' 



(22) 



respectively. As a result, /i^^-* will saturate to a balanced 
position at which the above two opposite currents cancel 
with each other. Because cos^ | — sin^ | = cos 6*, this 
balanced position of /i*-®^ turns out to be proportional to 
cos 6',^ 

The triangles in Fig. ^c) show the self-consistent re- 
sults of the NLSV signal as a function of 9. For compar- 
ison, the functions ±cos0 are also plotted by sohd and 
dashed lines, respectively. One easily concludes that the 
NLSV signal in the presence of 9 is given by 



.Rnlsv = ±^ cosf 



(23) 



where ± depends on the magnetization of the grain. 



C. Discussion 

There are several points should be clarified: 

(i) ^NLSV we obtained for the present device and pa- 
rameters happen to be of the same order as the experi- 
mental observations, where the NLSV signals are of the 
order of 1 - 10 mSlM. 

(ii) Above, we only consider the case that S decreases 
by 1/2 when the extra electron is added, where S is the 
magnitude of the total angular momentum of the grain. 
Also, according to Fig. [TTl there is small probability 
that S increases by 1/2 when adding the extra electron, 
i.e., the extra electron is preferred to be added to the 
majority band, then the spin selection rules will become 
totally reversed as 



(iVo,5'o,^.M,;t|iVo + l,5'o-f i + = 
(iVo,^o,5.Mu|A^o + l,^o + i = 



So + Sz + l 
25o + 1 '■ 
So-Sz + l 



25o + l 



(24) 



so that tunneling of spin-up (spin-down) electrons will be 
favored when S^ ~ 5*0 (S^ ~ —Sq). Therefore, the results 
will be totally reversed. This is a direct consequence of 
the strong Coulomb repulsion and the unequal spacings 
5a and 5i of majority and minority one-particle levels 
for an ultrasmall grain, and a major difference from the 
relatively large films i^i^ 

(iii) Note that in regions (4) and (5), the stead-state 
splitting of spin-up and -down chemical potentials is 
maintained by the spins continuously injected by Ic- This 
is different from in region (6), where there should be 
no current flowing in or out at steady state, because 
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of the voltmeter. Therefore, in regfon (6), the spin-up 
and -down electrons will eventually relax to one chem- 
ical potential for sufficient long time. That is why we 
consider, for the steady-state solution, only one spin- 
irrelevant chemical potential /i^^-* in region (6). 

(iv) In the simulation, although the parameters listed 
in Table U for a 5o = 1000 grain are exploited, we only 
use 5*0 — 100 for simulation because of the limited com- 
puting power. We have checked the results from So — 10 
through ^0 = 100, and the results turn out to be quanti- 
tatively unchanged as long as 5* 1. 



IV. MAGNETIZATION SWITCHING 

In this section, we will present the magnetization 
switching of the grain under a current Ic exceeding a crit- 
ical value. This critical Ic is determined by the gate volt- 
age Vg and the spin bias Vs , at which the strong Coulomb 
and magnetic blockades of the grain are lifted. Then cur- 
rents can flow through the grain, and transfer the angular 
momentum carried by the flowing electrons to the grain. 
This will be discussed in Sec. IIV Al 

Besides, still because of these blockades, one can not 
pre-assume there is always a pure spin current flowing 
through the grain, so whether the reversal is accompanied 
by the pure spin current is waited to be checked. This 
will be discussed in Sec. IIVBI 



A. Critical Ic and Vs 

As we have discussed, the electron tunneling between 
the lead and the grain are subjected to the Coulomb and 
magnetic blockades simultaneously. According to Eq. 
(|B12[) . by choosing suitable gate voltage, the charging 
energy Ec can be compensated, but the transition ener- 
gies still depend on the magnetization of the grain, i.e., 
the grain may be magnetic blockadediAS In this situation, 
this magnetic blockade can be lifted by applying the spin 
bias Vs exceeding a critical value, which thereby defines 
the minimal critical Ic and Vg. 

We will use the reversal from Sz — —5*0 to 5*0 to ex- 
tract the minimal critical Ic and Vg at which the mag- 
netic blockade can be lifted and the switching can be 
performed. 

Suppose the grain is initially prepared at the state 
\No, So, —So)- By adding a spin-up electron from lead 
(4,5) into the minority band of the grain, the grain will 
transit to the state \No + I, So — ^, —So + ^) ■ This tran- 
sition, energetically requires that 



(25) 



Via this transition, Sz increases by 1/2 unit. 

Further, by draining a spin-down electron from the mi- 
nority band of the grain to lead (4,5), the grain will 
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FIG. 7: (a) The spin bias Va, (b) the normalized magneti- 
zation, (c) the NLSV voltage, and (d) the NLSV signal as 
functions of the driving current Ic- The parameters 6 — 0, 
7 = 0.1, I^ = 0, So = 100, r = 20 mK, AVg = 0, and other 
parameters are the same as those in Fig. |3] 



transit from the state \No + I, So — ^,—So + 5) to 
jA'^Oj 'S'o, — + 1). This energetically requires that 



l,So- 



El 



Afo,So,-So-|-l)- 



(26) 



Via this transition, Sz also increases by 1/2 unit. 

Note that ^'"y^lL) = ±Vs/2. If one applied a sufficient 
large Ic, so that the spin bias Vg driven by Ic is large 
enough, for all the possible Sz £ [— 50,50], there are 
always 



(27) 

Then, one can expect a sequence of consecutive charging- 
discharging steps to be driven, which charges the grain 
with only spin-up electrons and discharges the grain with 
only spin-down electrons. As a result of this charging- 
discharging sequence, the magnetization of the grain will 
eventually be reversed from Sz = —So to Sq- 

Similarly, to reverse the magnetization from Sz = So 
to —So, the energy requirement is that for all the possible 

Sz 7 



(4), 



(28) 



Therefore 
^(^)(L)-..(^) 



the minimal required Vg, which equals 
fj,'"^' (L), is determined by the width of spectrum 
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El 



|JVo + l,So-i,S,±i> - ^\N o,So,S,) 

According to Eq. (|B12p . this spectrum width is 2Kn, 
which thereby set the value for the minimal required spin 
bias. The minimal critical current is then defined as the 
Ic by which the generated Vs > 2Kn. 

Fig. [7] shows Vs, magnetization, and NLSV signals as 
functions of Ic, when 9 — 0. For each Ic, Vs is calculated 
first using Eq. Then, Vs is put into the rate equations 
Eq. to self-consistently determine AV and n^^^ by 
using the same method shown in Figs. [5] (a3) and (b3) 
until both spin-up and -down currents through the grain 
vanish. Finally, the magnetization is obtained by putting 
the calculated /i'-^-' back to the rate equations. The NLSV 
signal i?NLSV is found by i?NLSV — AV/Ic- These steps 
is shown by Fig. O 

The triangles V and A in Fig. [7] indicate the P and AP 
cases we have already discussed in Fig. [SJ respectively. 
Keep increasing Ic until Vs exceeds 2Kn, the magnetiza- 
tion of the grain will be reversed. In the present set of 
parameters, the steady-state magnetization starts to re- 
verse when \Ic\ is a little larger than 2.5 mA. The switch- 
ing is accomplished after \Ic\ exceeds ~ 3 mA, at which 
\Vs\ is right larger than 2Kn = 0.02 meV. We attribute 
the broadening of reversal point at |/c| = 2.5 ~ 3 mA to 
the thermal fluctuation of the lead electron bath. At the 
reversal point, both the sign and slope of AV changes 
abruptly. i?NLSV also demonstrates a hysteresis loop in 
analogy to the hysteresis loop of the magnetization, but 
with opposite signs. This opposition has already been 
explained in the discussion (ii) of Sec. IIIII 

There are several points needed to be clarified: 

(i) The results should be qualitatively unaffected 
for small 6 ^ 0, because Q does not change 
-B|jVo+i,So-i,s.±i) - -E^IAfo.ScS,), while only these energy 
differences determine the critical spin bias Vs and Ic- 

(ii) We have concluded that the minimal critical cur- 
rent is only related to if at, which does not depend on 5*0, 
so we use 5*0 = 100 to perform the calculation. We have 
checked that the simulation results for other S'o ^ 1 turn 
out to be qualitatively unchanged. 



for all Sz € [— S'ojS'o]. the critical current density is about 



(ui) According to Eq. (|B9| . E\ 



E\ 



JVo + l:So-i,S,±i) 



JVo.So.Sj) is ^1^0 ^ function of the gate voltage Vg, 
therefore the critical Ic can be tuned by the gate volt- 
age. Above we only discuss the minimal critical Ic, i.e., 
the case when the Coulomb and the one-particle ener- 
gies are already compensated by choosing suitable gate 
voltage (details can be found in Appendix IB 5[) . There- 
fore, the spin bias only has to lift the magnetic blockade 
and thereby can be minimized. If the nearly degenerate 
situation were tuned away by a magnitude of AVg, the 
spin bias Vs then has to compensate the Coulomb and 
magnetic blockades simultaneously. Then, the critical Vs 
will become |2i4rjv| -f |Ayg|, which thereby requires larger 
critical Ic- 

(iv) The critical current density. For the present pa- 
rameters, the critical current is about /c = 3 mA, while 
the cross-section area is Aa = 170 x 65 nm^. Therefore, 



3 X 10-^A 



170 X lO-'^cm X 65 X lO-'^cni 



2.7 X lQ^k/(^m) 



This value is comparable to most experiments of nanopil- 
lars and multi-layers.— 



B. Pure spin current ? 



c 

O) 
O) 

o 



o 

"D 
CD 
> 

o 

03 
Q. 




time [1/r] 



FIG. 8: The spin-up and -down currents flowing through the 
grain and the magnetization of the grain as functions of time 
when assuming that the chemical potential of region (6) is: (a) 
/i'®' = Vs/2, (b) = -Vs/2- (c) region (6) is not connected 
to the grain. The positive direction of current is defined as 
flowing from lead (4,5) to lead (6). S = 0, K = S/Cat, 7 = 0.1, 
T = 20 mK, So — 100, AVg = 0, and other parameters are 
the same as those in Fig. [S] 



We are particularly interested in whether there is truly 
a pure current flowing through the grain during the re- 
versal in the present device. Therefore, we studied a 
situation that the grain is initialized at Sz = —Sq and 
Ic is suddenly switched on to generate a spin bias large 
enough to drive a switching from Sz — —Sq to 5*0, then 
see how the spin-resolved currents and the magnetization 
evolve with time. 

Roughly speaking, if we assume the electrons tunnel- 
ing through the grain transfer all their angular momenta 
to the grain, to reverse the grain from Sz = —Sq to 5*0, 
there should be at least 45*0 electrons tunneling through 
the grain during the reversal. Because the grain is much 
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smaller than region (6), it is safe to expect that the elec- 
trons flowing in and out region (6) along with the reversal 
process will hardly change /i'-^-' . 

According to Fig. [5] (a2), when Vg is positive while 
Sz = —So, ^^^^ will saturate at Vs/2; and when V^^ > 0, 
5"^ — So, /i^^^ will saturate at —Vs/2. So in the following 
we will compare two limits. The first limit assumes that 
;i(6) = Vs/2. The second assumes that = -Vs/2. 
These two limits are shown in Figs. [8] (a) and (b), re- 
spectively. 

Let us first consider the first limit /i^®-' — Vs/2 = 0.015 
meV. This situation is the same as shown in Fig. [5] 
(a2), but Vs is much larger in magnitude. According to 
Fig. m^a), as the magnetization changes from —Sq to 5*0, 
the spin-down current gradually becomes favored and the 
spin-up current unfavored, which is consistent with the 
spin selection rules Eq. p?| . Remember the configura- 
tion of chemical potentials remains unchanged during the 
reversal as shown in Fig. [5](a2). As a result, the magni- 
tude of spin-down current will keep growing during the 
reversal, and even after the magnetization is reversed to 
^0, there will still be a steady spin-down current flow- 
ing from lead (6) to lead (4,5). Although not shown in 
our result, one can expect, for sufficient long time, the 
spin-down current leaking from region (6) will eventually 
shift /i*^^^ down to —Vs/2, much like the same situation 
as shown by Fig. O (b2). Then the spin-down current 
will cease to flow. 

The second limit is shown by Fig. [5] (b). At the start, 
Sz = —So, so spin- up current is favored to flow from lead 
(4,5) to lead (6). As the reversal from —So to 5*0 goes on, 
the spin-up current will gradually become unfavored, and 
drop to zero after the reversal is accomplished. Although 
the spin-down current is favored when Sz = So, the Fermi 
levels for spin-down electrons on both sides of the grain 
are the same, so neither the spin-up nor -down current 
will be flowing after the reversal. 

After considering the above two limits, it is natural 
to expect that a real situation should be between them, 
i.e., /i*-^^ should float gradually from Vs/2 to —Vs/2, and 
there should be a spin-polarized current instead of the 
anticipated pure spin current flowing through the grain 
during and after the reversal, until the accumulation or 
drainage in floating region (6) is accomplished. 

The extra charge part of the current is produced along 
with the electron accumulation or drainage process in 
region (6). To prove this point, we decouple region (6) 
completely by letting F^g) — OM The results are shown in 
Fig. [5](c). In this case, the spin- up and down currents are 
always the same in magnitude and opposite in direction, 
i.e., there is only a pure spin current flowing. Besides, by 
integrating the current over time using the data of Fig. 

Qui = ilii) dt, (30) 

we obtain that both the spin-up electrons that tunnel 
into the grain and spin-down electrons off the grain are 



2Sq. This is consistent with the change of Sz by 2So 
during the reversal, where half is from the 2So incoming 
spin-up electrons and half the 25*0 outgoing spin-down 
electrons. 



C. Validity of approximation l'^ — 

Because surrounded by insulator, the grain is con- 
nected to the leads via the quantum tunneling. The ref- 
erence 20 of Ref. estimates that the tunneling rate 
for Ref. 8 is around F = 10^ s^"'^. According to Fig. [3 
where either the magnitude of spin-up or -down current 
is smaller than eF. Therefore, we estimate the current in 
and out the grain are well smaller than 

eF = 1.6 X lO^^^Coulomb x lO^s^^ = 1.6 x 10"^"v4(31) 

This value is 6 - 7 orders smaller than, e.g., the driven 
current Ic (~ 10^"^ A) shown in our Fig. [71 So it is valid 
for us to ignore in our numerical calculations, though 
we explicitly kept it in Eq. [51 



V. SUMMARY 

In this work, we theoretically studied the magnetiza- 
tion switching and detection of a ferromagnetic nanograin 
in a non-local spin valve (NLSV) device. 

Different from the original experiment, '^^ our 
nanograin is much smaller in size and at much 
lower temperatures, thus subjected to strong Coulomb 
and magnetic blockades. We describe the grain as a 
Stoner particle, whose ferromagnetism comes from the 
exchange interactions between itinerant electrons inside 
it. Because of the ultrasmall size, the one-particle levels 
inside it are quantized. Because of the ferromagnetism, 
the level spacings for the majority and minority electrons 
are unequal. 

As shown in Fig. [H^a), the grain is coupled to regions 
(4)-(6) of the NLSV device via quantum tunneling. Re- 
gions (4) and (5), and region (6) can be regarded as two 
nonmagnetic leads, respectively. In the lead formed by 
regions (4) and (5), a spin-dependent splitting of chemi- 
cal potentials (spin bias) is induced by the spin-polarized 
current Ic injected from ferromagnetic injector (1), as 
shown in Figs. [Sfa) and [3Kb). The other electrode (6) is 
floating and connected to region (5) through a voltmeter. 

By applying a Ic and measure the voltage difference 
Ay between regions (5) and (6), the magnetization of 
the grain can be read out by the NLSV signal i?NLSV = 
AV/Ic. Because of the unequal level spacings for the ma- 
jority and minority electrons and Coulomb blockade, the 
NLSV signal in the tunneling regime depends not only 
on the magnetization of the grain, but also on whether 
the majority or minority band of the grain is favored to 
contribute to the electron transport. The results when 
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the minority band is favored are right opposite to when 
the majority band is favored. In the presence of an an- 
gle between the easy-axis of the grain and the spin- 
quantization direction of the electrode, the NLSV signal 
is proportional to cos 6* and vanishes when 6 = 7r/2. 

By applying Ic exceeding a critical value, the magneti- 
zation of the grain can be switched reversibly by the spin 
bias generated by the Ic- Because of the strong Coulomb 
and magnetic blockades, the electron flowing between the 
grain and the electrodes is not possible unless both block- 
ades are lifted, then the angular momenta carried by the 
flowing electrons can be transferred to the grain. There- 
fore, the critical value of Ic to drive the magnetization 
switching is determined by: at what gate voltage Vg and 
spin bias Vg, both the Coulomb and magnetic blockades 
can be lifted. We also show that the current accompa- 
nying the switching may not be a pure spin current, due 
to the accumulation or drainage of electrons in the float- 
ing lead used for the NLSV measurement. A possible 
solution is to remove the floating lead. 

Our numerical evaluations using realistic parame- 
ters from the recent NLSV^i^i and the cobalt grain 



experiment* 



.9.10.11 



show that it is possible to employ the 



NLSV device to detect and switch the magnetization of a 
ferromagnetic nanograin under the present experimental 
conditions. 
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APPENDIX A: BOUNDARY CONDITION (III) 
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gram 



With the help of these arrows, one can write 

+ (Al) 

i.e., the current flowing into the node equals to those flow- 
ing out. The current density is related to electrochemical 
potential by Jo- = —(s'a I &)dxPta- We will take the spin-up 
component as a example. In (4) and (5), CTj = and 
the spin-up current flowing from the node to the grain is 
defined as 1^ = 1'^. Put these together, one can easily 
obtain that 



Ao^d.^^\'\L) = I^ 



"2^ 



(A2) 



Similarly, for spin-down /J^ — —I^, so 



'2^ 



A, 



'2^ 



^x^if{Q) (A3) 



APPENDIX B: THE QUANTUM THEORY OF 
THE FERROMAGNETIC NANOGRAIN 



1. Model of ferromagnetic nanograin 

In this work, we will describe the ferromagnetic 
nanograin by using the minimal possible model^i^ 
proposed to describe the experiment transport spec- 
tra through cobalt nanograin. ^'^ This model was then 
adopted to discuss spin-polarized current induced relax- 
ation and spin torque f^^iiii^ It also provided a starting 
point to study the Kondo resonance in the STM spec- 
trum of a ferromagnetic cluster on metal surface.— A 
more detailed microscopic tight-binding model with ex- 
change interactions and atomic spin-orbit couplings was 
also proposedfSi^i^ to reveal a unified origin of the mag- 
netic anisotropy as well as collective and quasiparticle 
excitations in the ferromagnetic nanograin. Besides, the 
Jaynes-Cummings model also reproduced the transport 
features by considering the interaction between particle- 
hole excitation and magnon. — 

In this work, we focus on how the collective spin and 
one-particle excitations of the grain react to the spin bias, 
therefore the minimal model is adequate for the current 
topic. The Hamiltonian for the grain and its couplings 
to nearby leads is given by 



total 



G 



(Bl) 



FIG. 9: The zoom-in of Fig. [TJb) near interface (4)/ (5) /grain. 
Arrows mark the locally defined positive direction in each 
region. Three regions meet at the central node. In the local 
coordinates of (4), the node is at L. In the local coordinates 
of (5), the node is at 0. 

The boundary condition can be found with the help 
of Fig. [HI which is a zoom-in of Fig. 1(b) near inter- 
face (4)/(5)/grain. The positive direction of the locally 
defined coordinate in each section is marked by arrow. 



where the Hamiltonian for the ferromagnetic grain takes 
the form, 

= V e,Aad^. - ^ S • S - + EcSN^ + VgSN 



(B2) 



where the first term stands for the kinetic energy of elec- 
trons in the grain, dir, {d}^^ ) annihilates (creates) an elec- 
tron on the one-particle level i in the grain, with energy 
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e^a and spin a = {T, i}. S = ^ I]<t<t' Aa^ <y<y"^i'y' is 
the total angular momentum of the grain electrons, 't 
is the vector of Pauli matrices. J is a phenomenologi- 
cal constant depicting the exchange interactions between 
each pair of electrons in the grain. Sz is the z-component 
of S. if AT is the volume-independent anisotropic con- 
stant. In this work, we consider that the fluctuation of 
the electron number 1) is much smaller than than the 
itinerant electrons > 1000 already in the grain, so the 
fluctuation of Kyi as a function of the electron number is 
ignored. bN is the number of extra electrons added into 
the grain, compared with a reference electron number A^o 
already in the grain. Eq is the charging energy required 
to add the excess electrons in the grain, which we will see 
later can be compensated by applying a gate voltage Vg. 

We refer the part where regions (4), (5), and (6) con- 
nect the grain as two "leads" , one is from regions (4) and 
(5) together, and the other from region (6). For conve- 
nience, we call them lead (4,5) and lead (6), respectively. 
The Hamiltonian for the leads takes the form 



/c,a,T 



(B3) 



where cl.^^{ckar) is the creation (annihilation) operator 
for a continuous state in lead a € {(4, 5), (6)} with energy 
eka and spin t G {-|-, — }. The tunneling between the 
grain and the leads is described by 



k,a,i 



+ (sin -fcl^+ + cos "-fcl^Jd.i] + H. 



c, 
(B4) 



where 9a is the angle between the easy axis of the grain 
and the spin-quantization direction of lead a. For sim- 
plicity, we set the easy-axis of the grain as z-axis and 
assume that ^(g) = and ^^(4^5) = 9 £ [0,7r/2], as shown 

in Fig. m In the following, we denote = when 

^(4,5) = for simplicity. 



2. Ground branch \No, So, Sz) of the grain and 
possible excitations 

The eigen states of Hq are labeled by \{ni},S,Sz), 
where {ui} denote the occupation on each level in the 
grain, S and 5*2 are, respectively, the quantum numbers 
for the magnitude of S and Sz- Because of the low ex- 
perimental temperature (as low as 20 niK) and the large 
Coulomb repulsion Ec (> 30 meV),^ the charge fluctu- 
ation and particle-hole excitations [Fig. (c)] are sup- 
pressed hence it is reasonable to assume that the elec- 
trons in the grain compactly occupy all the lowest avail- 
able levels. These states are denoted as \N,S,Sz), where 
A'^ = Ui is the total electron number in the grain. 



(a) 



-©- 



-e- 



majority 



-e- 



minority 



(b) 



8N 



majority 



8N. 



minority 



(c) particle-hole 
^9 



I! 



(d) spin-accumulation 




(e) one-particle 

— 

-Q 



-e- 



FIG. 10: The filling of one-particle levels in the majority and 
minority bands of the grain with (a) A'^o electrons, and (b) 
iVo -I- SN electrons, where 5N = 5Na + 5Ni, and 5Na and SNi 
are the excess electrons added to the majority and minority 
band, respectively. For the ground branch with A^o electrons 
in the grain, the magnitude of the angular momentum S — Sq. 



Ep°^, and Epl are the highest occupied levels of the majority 
and minority bands for the ground branch for N = No and 
S = So- Sa and Si is the level spacings near the and 
Ep"-, respectively, (c) The particle-hole excitations, with ex- 
citation energy 5i,5a > ImeV. (d) The spin-accumulation 
excitations, with excitation energy ~ J ~ 2meV. These two 
excitations are high-energy excitations, with the excitation 
energy larger than the spin bias Vs, thus can not be excited 
and will be omitted. Note that according to Fig. |3] Vs ~ 0.1 
meV. (e) One-particle excitations that changes the electron 
number by one. It can add extra electrons to the majority 
or minority band. By tuning the gate voltage, its excitation 
energy can be minimized down to ~ Kn ~ O.OlmeV, which 
is smaller than the spin bias Vs ~ O.lmeV. Then, one-particle 
excitations is the only possible excitation for the present de- 
vice and parameters. 



When there is A^ = A'o electrons in the grain, the com- 
petition between the kinetic energy and the J term will 
force No/2 + So electrons (the majority band) to ori- 
ent anti-parallel with the rest Ao/2 — 5*0 electrons (the 
minority band), leading to a nonzero magnitude of an- 
gular momentum S" = 5*0 at the ground states(Stoner 
instability), as shown in Figs. S] and [TO] (a). This 
overall ground branch is denoted as \Nq, So, Sz) , where 
Sz e [-S'o,5'o]. The {2So + l)-fold degeneracy of the 
overall ground branch is lifted by the anisotropy, with 
two degenerate ground states \No, So, ±S'o)- 

There are three kinds of basic excitations from the 
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ground branches \No, So, Sz), as shown in Fig. [roTcl. (d), 
and (e). The particle-hole excitation destroys the com- 
pact occupation of the one-particle levels, while does not 
change N, S, and Sz. The spin-accumulation excitation 
changes S by moving electron between the majority and 
minority bands, while does not change A'o and Sz- The 
one-particle excitation changes N by adding or removing 
electrons, which also leads to changes in S and Sz- 

Generally, the excited energy from |A^0) <S'o, '§'2) to | Afo + 
SN, So + 6S, Sz + SSz) can be found with the help of Figs. 
[TOTa) and (b) as (using the relations SN = 6Na + SNi and 
SS = {SNa - SN,)/2)^ 



5E{5N, SS, SSz) 

El 



No+5N,So+SS,S^+SS^) 

{SNf[Ec + + Is^] + SNiEp 
o o 



1. 1. T.l 

jSa + -6^ + Vg] 



(SS) 



6NSS[-Sa 



2 



SSi^ 



2 



K 



N 



{2Sz5Sz + SSl), 



(B5) 



where the definitions of parameters are given in Fig. [10] 
and eI' ^{E'^\ + eI"^/2- 



The stability of the branch \N = No, S = So) requires 
that its energy should be at least locally minimal com- 
pared to the states \N — No, 5 = 5o ± 1), i.e.. 



SE{SN ^0,SS = +1) = Sa + - J{So + 1) > 0, 
5E{5N = Q,5S= -1) = Si - + JSq > 0, (B6) 



where we have used Eq. (jBSP and ignored the anisotropy 
term because Kn is much smaller than other energies. 
In other words, Eq. (|B6p leads to that, for J belongs to 
the range 



A So 



So 



< J < ^ 



So 



(B7) 



the grain will adopt a stable S = So when there are Nq 
electrons in the grain. Once J exceeds this range, the 
overall ground branch will evolve to adopt a smaller or 
larger S* = = 5o T 1, then Eq. dBT]) still holds for the 
new S'o (note that A'p"^^ = Af ±6aTS^). 



One-particle excitations to the branches 

|iVo-fl, So ±1/2,5'.) 



TABLE 

energies are in meV 
~ '.10.11 



I: Parameters for 

45.61 



Sa, 



a So 
Ec, So 
S 



by experiments,- 
calculations.— According to Eq, 
into Vg, thus is set to for convenience. ^' 
zero point throughout the paper 



-- 1000 grain. All 
Kn are estimated 
are from the band 

Ep can be absorbed 



(5) 



(cx)) is the energy 



Ec 


Kn 


5a 


Si 


A So 


Ep 


^t(5'(c») 


30 


0.01 


4.61 


1.19 


2000 









We will adopt a set of parameters for a grain with 
5*0 = 1000, as given in Tab. HI The theoretical 
calculations"*^'^"^ based on these parameters are con- 
sistent with most features observed in the experimental 
transport spectrai^ii^ With these parameters, the value of 
J can be deduced using a saddle point analysis as follows. 



Range of ,7 for a stable 5a by saddle-point 
analysis 



(A^"-&)/S„ (A^"+5^-&)/(S„+1/2) (A^"+5j/(S„+1) 



■§ 1001.0- 

E 1000.5- 
o 

E 1000.0- 

? ! 

15 998.5^ 



(a) 

1 


No-electron | 






1 

(No-^1)-electron 


f 

______ J 






1.996 



1.998 



2.000 

J [meV] 



2.002 



2.004 



2.006 



Above we just assume that S = So when there are N = 
No electrons in the grain. Since So originates from the 
competition between the J term and the kinetic energy, 
its value should be calculated for the given values of J and 
5a,i- However, since we already know Sq 1000 from the 
experimentsr"^ i^°'^-^ and 5a,i from the band calculatiouf^ 
we will use these data to conversely deduce the value 
of J, then generate other information (e.g., one-particle 
excitations) at the proximity of the deduced J, much like 
a saddle-point analysis. 



FIG. 11: (a) Total angular momenta for the ground branches 
of the grain with A^o and A'o -I- 1 electrons, respectively, (b) 
Along the zigzag line of (J, Vg), the A'o- and (A'^o -I- l)-electron 
ground branches are nearly degenerate. The parameters are 
listed in Table HI 

The excitation energies for the particle-hole excita- 
tion shown in Fig. [TU] (c) and spin-accumulation exci- 
tation shown in Fig. [in] (d) are of the order of Sa/i and 
J, respectively ! ^^'"^^1^^ According to the parameters for 
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So = 1000 shown in Tab. HI these excitations are of meV, 
much larger than the Vg ^ 0.05 meV estimated in Fig. [31 
thus wiU be excluded. In the following, we will only con- 
sider the one-particle excitation from the A'q to iVo -I- 1 
electrons, as shown by Fig. [10] (e), i.e., adding excess 
electrons to the grain. 

Again, because Ec ^ ksT, the gate voltage in this 
work is restricted so that one and only one excess electron 
{SN = 1) can be added into the grain. Depending on this 
excess electron being added to the majority or minority 
band, the magnitude of the angular momentum could 
change to S' = S'o + 1/2 or — 1/2, respectively. With 
the help of Eqs. (|B5|) and (|B7|) . we find that, when 



aSq 



Sn 



< J < — 



(B8) 



^0 + 1/2 ' 

the iVo -t- 1-electron ground branch will adopt S — 
Sq — 1/2, and the required transition energies from 
\No,So,S,) are 



No+l,So- 



V,+E, 



El 



\No..So..S,) 



C 



fpSo 



J 



{So 



(B9) 



On the contrary, when 

6a - 5i 



A So 



< J < 



(BIO) 



^0 + 1/2 - - So + 1 

the A^o + 1-electron ground branch will adopt S ~ 
So + 1/2, and the required transition energies from 
\No,So,Sz) are 



E, 



El 



J 



Kn,1 



(Bll) 

To summarize the above saddle-point analysis, the 
angular momenta of iVo- smd [Nq -j- l)-electron ground 
branches as a function of J are shown in Fig. [TlT a). 
in which the arrows mark the ranges indicated by Eqs. 
(iBTl) . (iB8l) and (iBlOl) . 

According to Fig. fTTT a). for most value of J, S of 
the {No + l)-electron branch is 1/2 smaller than that 
of the A^o-electron branch. This is a direct results of 
6i < 6a- Although J is not a tunable quantity, Fig. 
Illf a) implies that in reality the excess electron is far 
more likely to occupy the minority band than the ma- 
jority band, which is also consistent with the previous 
literatures. ^^'^ Therefore, in the following we will mainly 
consider the one-particle excitations between the ground 
branches |A^o, So, S^) and |7Vo 1, - 1/2, S:,). 



5. States used for numerical simulations 

Remember we have set ii^^\oo) as the energy zero 
point. With respect to fi^^\oo), we can always choose 



suitable Vg in Eqs. (|B9p and (jBll[) to compensate Ec 
and other energies, so that the ground-branches with 
A^o and A^o + 1 electrons can be tuned to be nearly de- 
generate. These as a function of J are shown in 
Fig- fllTb). For instance, in Eq. (|B9p . by choosing 
K - - 



Er 



obtains 



^Fi 



J{So + l/4)/2 -I- A'Ar/4S'o, one 



El 



|Aro + l,So-i,S,±i> 



El 



\No,So,S,) 



T^S.. (B12) 



In this context, the energy required to add an electron 
from lead (4,5) into the grain is related to only the mag- 
netization of the grain. According to Eq. (|B12[) . the spec- 
trum width of El 



-E\No,So,S:,} is 2Kn for 



all the possible Sz G [—So, So]- This value is the key to 
determine the critical Ic in the Sec. IIV Al 

Experimentally, it is easy to find the suitable Vg at 
which the ground-branches with No and Aq -|- 1 elec- 
trons are nearly-degenerate, as in the usual transport 
experiments For example, one just apply a small 
charge bias voltage > 2Kn between lead (4,5) and lead 
(6), and measure the current through the grain while 
scanning Vg, like a usual source-gate-drain measurement. 
Because of the Coulomb blockade, the grain can not 
conduct electrons unless the A'o- and (A'o -I- l)-electron 
ground-branches are degenerate. Therefore, the nearly- 
degenerate situation is find as: at which Vg, the grain 
is conducting under a small charge bias voltage between 
(4) and (6). 

Based on the above analysis and discussions from 
Appendixes IB 2IIB 5[ in the following, we will consider 
mainly the one-particle excitations from the branches 
lA^o, So, Sz) to |A^o + l,So- 1/2, Sz) and when then en- 
ergies are nearly degenerate. Specifically, we will choose 
a set of {J,Vg) marked by the circle in Fig. [TTT b) for our 
numerical simulations, where J = Ap° /{So + 1/2), and 



Vg = -S, -Ec- Epl - J{So + l/4)/2 + Kn/^So + AVg, 
AVg is a small variation of the gate voltage that drives 
the grain away from the nearly-degenerate point of A^q- 
and {No + l)-electron ground branches. 



APPENDIX C: VALIDITY OF RATE EQUATIONS 

Although the rate equations formalism is widely em- 
ployed for the mesoscopic systems weakly coupled to the 
electrodes, its validity deserves some discussioni^^i^ 

In the previous works by Waintal et al^^^ the in- 
trinsic spin relaxation is considered in terms of coupling 
to a bosonic bath. We do not consider this effect for 
two reasons: (i) According to Eq. (3) of Ref. [l^, the 
intrinsic relaxation will lead to a term similar to the 
Gilbert damping, which tends to relax the grain to one of 
its two degenerate maximally magnetized ground states, 
e.g., \No,So,-So) or \No,So,So). In the following, we 
will mainly discuss how to use the NLSV signal to read 
out the magnetization of the grain (Sec. IIIip . which al- 
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ready limits the discussion to these maximally magne- 
tized states. Therefore, the results in Sec. IIIII will be 
qualitatively unaffected by the intrinsic relaxation, (ii) 
For the spin bias -induced magnetization switching dis- 
cussed in Sec. IIVI the results are only valid when the 
coupling between the grain and the lead is dominant and 
much smaller in time scale than the intrinsic relaxation. 
On the other hand, the Born and Markoff approxima- 
tions enforce that T is much smaller than ksT (bath 



temperature) and the energy difference between many- 
body states in the grain. These two requirements confine 
the validity range of F in this work. 

Besides, the particle-hole excitation terms in Eq. (10) 
of Ref . are also absent in this work, because the energy 
scale of these excitations are of the order 6a, Si > ImeV, 
which already one order larger than the spin bias {Vg ~ 
O.lmeV) that can be generated according to Fig. [S] 
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